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0^ : ABSTRACT 

^ ■ HST observations of some QSOs show a strong, abrupt increase in polarization at 

. rest wavelength about 750 A. The closeness of the polarization rise to the H I Lyman 

psj , edge suggests a connection, but the displacement to shorter wavelengths, and the shape 

I of the polarization rise require explanation. 

I We have computed the polarized spectrum of a thermally emitting accretion disk 

I around a supermassive black hole, including the effects of the relativistic transfer func- 

i-C ■ tion. The local stellar atmosphere spectra show a blueshifted polarization rise in the 

Lyman continuum, as found by Blaes and Agol (1996). However, the relativistic trans- 

2 ■ fer function adds an additional blueshift of sufficient magnitude that the model cannot 

-4-^ I 

c/3 . explain the observations. 

I We show that a good fit results if the emitted radiation is assumed to have a sharp 

I increase in polarized flux at the Lyman edge in the rest frame of the orbiting gas. 

^ I Relativistic effects then cause the observed polarization to rise sharply at a wavelength 

' substantially less than 912 A. The blueshift depends on the angular momentum of 

the black hole and the inclination of the disk. A good fit to PG 1630+377 results 
from a simple model with a dimensionless angular momentum a* = cJ/M^ = 0.5 and 
an observer viewing angle /Uq = cos Oq = 0.1. Smaller values of a* give insufficient 
blueshifts, and values close to = 0.9982 require unrealistically large polarizations in 
the rest frame of the gas. An intermediate value of a* might result from coallescing black 
holes, successive accretion events, or electromagnetic extraction of angular momentum 
from the hole. 



Subject headings: galaxies: active — quasars: general — accretion, accretion disks — 
polarization — black hole physics 
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1. Introduction 

The leading model of energy production in QSOs involves an accretion disk orbiting a super- 
massive black hole (Blandford 1990). Evidence for black holes in the nuclei of nearby galaxies is 
increasing (Rees 1997), but proof of accretion disks in QSOs has been elusive. 

The expected disk effective temperatures imply disk emission at wavelengths in rougli agree- 
ment with the Big Blue Bump observed in energy distributions of QSOs (Shields 1978; Malkan 
1983; Czerny and Elvis 1987). Accretion disk fits to the energy distributions of a number of QSOs 
have been quite successful (e.g., Sun and Malkan 1989, Webb ct al. 1993). These models indicate 
masses Mjj of over 10^ Mq and accretion rates M above 1 yr^^ for luminous QSOs. Lyman 
edges in the total flux should be inconspicuous because of relativistic Doppler broadening (Laor and 
Netzer 1989, hereinafter LN; Laor 1991) and NLTE effects (Coleman 1993; Shields and Coleman 
1994, hereinafter SC; Hubeny and Hubeny 1997). The Lyman edge of H I will be most conspicuous 
in relatively cool disks {Tmax — 25, 000 K), but will appear as a slope change over the wavelength 
range ~ 1100 to 700 A (SC) rather than an abrupt discontinuity. 

Models of accretion disks predict substantial electron scattering opacity at optical and ultravi- 
olet wavelengths, leading to polarization of the emitted light (see section 2.1). In contrast, observed 
polarizations in the optical and near ultraviolet arc typically ~ 1 percent or less (Stockman, Moore, 
and Angel 1984; Webb et al. 1993). Moreover, in radio loud objects where a radio extension in- 
dicates the position angle of the disk axis, the observed polarization tends to be parallel to the 
disk axis (Berriman et al. 1990). Electron scattering in a geometrically thin disk would produce 
polarization perpendicular to the axis. Predicted polarizations can be reduced by absorption opac- 
ity, Faraday rotation (Agol and Blaes 1996), general relativistic depolarization (Laor, Netzer, and 
Piran 1990, hereinafter LNP), and disk surface structure (Coleman and Shields 1991); but some 
concern remains (Antonucci et al. 1996). 

LNP computed the expected energy distribution and polarization of QSO disks with a simple 
model of the atmosphere and full consideration of the general relativistic transfer function. They 
found that while the Lyman edge in the energy distribution is inconspicuous in general, a substantial 
drop in polarization in the H I Lyman continuum should result from the increase in absorption 
opacity relative to scattering. Several groups undertook spectropolarimetry with the Hubble Space 
Telescope (HST) Faint Object Spectrograph (FOS), targeting "Lyman edge candidates" for which 
International Ultraviolet Explorer (lUE) spectra suggested an intrinsic drop in the flux around the 
Lyman edge (e.g., Koratkar, Kinney, and Bohlin 1992). Rather than a drop in polarization in the 
Lyman continuum, these observations revealed an abrupt rise in polarization at rest wavelength 
~ 750 A in several objects (Impey et al. 1995; Koratkar et al. 1995). The polarization increases 
by a large factor, from < 1 percent at A > 1000 A to ~ 5 percent in PG 1222-1-228 at A ~ 600 A 
and ~ 20 percent in PG 1630-1-377. The proximity of the polarization rise to the Lyman edge is 
suggestive (Impey et al. 1995). 

An intriguing explanation of the displaced polarization rise was offered by Blaes and Agol 
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(1996, hereinafter BA). Using LTE and NLTE models for the disk atmosphere with a hydrostatic 

vertical structure, they found that, for Tejj ^ 20,000 K, a similarly blueshifted, steep polarization 
rise resulted from a combination of effects. The polarization drops at the Lyman edge because of the 
increased absorption opacity, but then it increases with decreasing wavelength because of the 
dependence of the H I photoionization opacity and the increasingly steep source function gradient 
at higher frequencies. A strong source function gradient leads to severe limb darkening, needed for 
strong polarization (Agol, Blaes, and lonescu-Zanetti 1997; Cheng et al. 1988). A model for a disk 
around a nonrotating black hole gave a promising fit to PG 1222+228, but polarization as strong 
as observed for PG1630+377 could not be achieved. However, BA did not include the relativistic 
transfer function, and it was unclear whether these effects would upset the fit to PG1222+228. 

We report here results of a theoretical study of polarization in QSO accretion disks. We 
focus on the steep polarization rises observed in PG 1222+228 and PG 1630+377. In Section 2, 
we describe our stellar atmosphere models and the treatment of the relativistic transfer function. 
Then we describe attempts to fit the two QSOs with accretion disk atmospheres, and show that 
general relativistic effects preclude an acceptable fit. In Section 3, we show that a good fit can be 
achieved in a simple model with a sharp polarization jump at the Lyman edge in the rest frame 
of the gas, and we consider possible sources for this emission. Our conclusions are discussed in 
Section 4. 



2. Accretion Disk Model Atmospheres 

2.1. Computational Methods 

We have calculated LTE stellar atmospheres for conditions appropriate for QSO disks. These 
models were supplemented by some "mixed" LTE-NLTE calculations as described below. Because 
we arc interested only in the continuum, and because a perfect description of the vertical structure 
of AGN disks is not available, we used simple methods to model the atmosphere. These resemble the 
LTE calculations of Coleman (1993), but the codes used here were written largely independently. 
The physics included is similar to the treatment of BA, but the computational methods are quite 
different. 

Atmospheres were computed for a set of radii and corresponding values of effective temperature, 
Tgff , and vertical acceleration of gravity, gz , evaluated for a given Mh, M, and angular momentum, 
a* = cJ/GM"^, from the relativistic relations given by Novikov and Thorne (1973) and Page and 
Thorne (1974). The atmosphere was computed for a set of points (typically 91) equally spaced 
in the logarithm of the column density (gm cm~^) ranging from log S(z) = —6.5 to 2.5. The 
hydrostatic density profile was integrated from the top down according to the equation 



dPg/dL = 5e// =gz- grad, 
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where Pg is the gas pressure and grad = npF/c is the radiative acceleration. Here 

/•oo roo 

Jo Jo 

is the flux- weighted mean opacity per unit mass (Mihalas 1978). 

Radiative transfer in the atmosphere was calculated by performing a downward integration of 
the equations 

dKy/dxy = Hj, 

and 

dHy/dTu = €y{Jy - By), 

where , H,j , and K^y arc the usual Eddington moments of the radiation field Iy{Ty,^), Ty is 
the monochromatic optical depth from the top of the atmosphere, = K^/^Ky + a^), and and 
ay are the absorption and scattering opacities per unit mass. The emergent flux Hy was adjusted 
in a shooting method to satisfy a lower boundary condition J^, = at a, fiducial optical depth, 
typically = 8. The first iteration used a gray atmosphere temperature profile, the Eddington 
approximation (J^ = SK^ ), and a surface boundary condition = \/?,Hy. Subsequent iterations 
used variable Eddington factors and a temperature profile corrected by the Unsold-Lucy method 
to achieve flux conservation (Mihalas 1978). We assumed constant flux, i.e., no local heating (see 
below). The resulting {ti, ) was used to compute the source function Si, {td ), and new values 
of Ju , Ku , and were then computed at each depth point from the usual exponential integral 
expressions (Mihalas 1978). These gave new values for the Eddington factors, the quantities needed 
for the temperature correction, and Kp . The entire procedure (hydrostatic integration, transfer 
integration, moments, temperature correction) was iterated as necessary, typically three iterations. 
Opacity sources included were electron scattering, free-free absorption, and bound-free absorption 
from the lowest 7 levels of H I and the ground states of He I and He II. Ionization and H I 
level populations were computed from the Boltzmann-Saha equation (LTE). Convergence was good 
except in a thin surface layer for certain ranges of T^ff , for which H or He were partially ionized 
at the surface. The flux in the related Lyman continuum of H or He was negligibly small in these 
cases, and the fluxes at longer wavelengths were not significantly affected. 

The assumption of constant flux through the atmosphere should be adequate in our case (see 
also BA). We describe below models using the same parameters (Mij,M, a) as BA. The region of 
the disk emitting the Lyman continuum is characterized by the radius 9.5GM/c^ giving maximum 
Tf,ff (25,500 K). In the vertically averaged model of NT with a = 0.1, the surface density of the 
disk (midplane to surface) is He = 2000 g cm~^. The surface density above the point where 
= 1 at A = 911 A is Si = 0.8 g cm . If energy is dissipated along the vertical extent of the 
disk roughly in proportional to density, then only a small fraction of the flux is dissipated in the 
photosphere. Also, we find that, at Ti, = 1, the dissipation rate per unit mass is tiny compared with 
the radiative heating and cooling rates in the atmosphere. These conclusions are even stronger 
if dissipation is proportional to gas pressure or total pressure. Moreover, our conclusions depend 
more on relativisitic effects than atmospheric details. 
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Polarization was calculated using a straightforward iterative procedure by Voigt (1951) applied 
to the final converged atmosphere. Although Voigt studied cases with small scattering opacity 
fractions, we found that convergence was rapid in all cases. The exact polarization exceeds that 
of the first iteration by almost a factor two for wavelengths at which electron scattering strongly 
dominated the opacity. This results from multiply scattered photons that produce substantial 
"prepolarization" of the radiation going into the last scattering. 

CS and others have emphasized that the effective gravity ge/f in AGN disks will be very low 
because of close cancellation of gz and grad ■ This leads to a large gas pressure scaleheight, a 
low density in the photosphere, high ionization, and a large opacity contribution from electron 
scattering. Coleman (1993), Hubeny and Hubeny (1997), and others have shown that significant 
differences in emitted spectrum can result from use of a fully self-consistent vertical structure rather 
than an approximate, constant gravity solution. However, our results show that, for the polarization 
phenomena of primary concern here, the details of the local atmospheric emission are less important 
than the effect of the relativistic transfer function. Therefore, we followed the procedure of BA 
in using a value of gz that was 20 percent larger than the Eddington limit for a given T^ff . 
This is about the minimum gz that allowed convergence for most T^jf . For these purposes, we 
took gEdd = '^es'^F/c, where k^s is the electron scattering opacity for fully ionized pure hydrogen. 
(Reduction of gz to l.lgEdd, in cases where the atmosphere still converged, gave considerably lower 
values of geff , but the emergent spectrum and polarization changed only modestly.) 

Sensitivity of our results to NLTE effects was evaluated in several sample cases by using our 
LTE atmospheric structure, p{T,) and r(S), as input to an NLTE stellar atmosphere program 
developed by the Munich University Observatory stellar atmosphere group. This code, using mod- 
ern methods including accelerated lambda iteration (ALI), includes bound-bound and bound-free 
transitions of many levels of H and He. When set to run in LTE mode, this code gave emergent 
spectra in excellent agreement with our LTE code. In NLTE mode, the results showed the expected 
reduction of the Lyman absorption edge for relatively cool atmospheres (cf. Coleman 1993; Hubeny 
and Hubeny 1997). However, the polarization as a function of wavelength was not greatly different 
between the LTE and NLTE cases, in agreement with BA. 

Following BA, we cut off the disk at 50Rg, where Rg = GM/(? is the gravitational radius. For 
Schwarzschild black holes, we computed atmospheres at r* = R/Rg = 7.76, 9.52, 15.3, 20.3, 31, 
44, and 50. Test cases with additional radii gave similar results. Summation over disk radii and 
application of the relativistic transfer function, which accounts for Doppler shifts, transverse and 
gravitational redshifts, and gravitational deflection of light rays was accomplished with a modified 
version of the computer codes used by LNP and kindly made available by T. Piran and A. Laor. 
These codes use a set of 51 radii from r* = 1.2346 to 400 to cover disks around either rotating 
or nonrotating black holes. We typically used as set of 80 frequencies uniformly spaced in log v 
between 1.96 x 10^^ to 1.85 x 10^^ Hz. 

The LNP transfer function assumes limb darkening laws in the intensity and polarization for 
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pure electron scattering (Chandrasekhar 1960), although there is provision for a scale factor in the 

polarization. The actual limb darkening in the polarized intensity can be quite different (cf. Cheng 
et al. 1988; Agol, Blaes, and lonescu-Zanetti 1997; and references therein). In order to accomodate 
a general limb darkening law with minimum changes to the LNP codes, we used a simple fitting 
procedure. The polarized intensity from a given radius and frequency can be described in terms 
of the orthogonal polarized intensity components ir-(/^) and /;(/x). Each of these was fit with a six 
term series in integer powers of 1 — from n = to 5, which gave an excellent fit in all cases. 
Heuristically, this may be thought of as a set of 12 accretion disks, each of which emits at its 
surface either pure Ii or pure 7^, with a limb darkening (1 — ^)". The corresponding 12 transfer 
functions were computed with the LNP code and applied to the respective fitting components, and 
the resulting Stokes parameters at the observer were summed to produce the observed spectrum. 
Our atmosphere results were interpolated to the LNP radius grid and the adopted frequency grid 
using a linear interpolation in R and v for (1) the brightness temperature defined by -Fj^ = B^{T^) 
and (2) / Ir{^) and Ii{p) / li{0) at five values of n (prior to determining the fitting coefficients). 
Test cases gave good agreement with a disk emitting with the electron scattering darkening law 
and with results from an independent transfer function code by Agol( 1997a, b). 

The emergent energy distribution and polarization for an atmosphere with Tgjj = 25, 500 and 
log gz = ^og{1.2gEdd) = 2.58 are shown in Figure 1. The polarization rises in the Lyman edges 
of H I and He I are evident, similar in the case of H I to the results of BA (who did not include 
He). Note also the substantial polarization at optical wavelengths (see also Figure 3 of BA). The 
low effective gravity leads to a low density at the optical photosphere, and electron scattering 
contributes substantially to the opacity ('^ 75 percent at Ty = \ for A = 5000 A). Some earlier 
models (e.g., LNP) used the midplane density throughout the atmosphere. The predicted optical 
polarizations may be a problem for disk models, given the weak observed optical polarizations for 
radio quiet QSOs (e.g., Berriman et al. 1990). 



2.2. Results 

We have computed models for PG 1222+228 and PG 1630+377, for which the most dramatic 
Lyman continuum polarization rises have been reported. For PG 1222+228, following BA and 
Webb et al. (1993), we have used a nonrotating black hole with Mg = Afj^/lO^ Mq = 5.3 and Mq = 
M/1 Mq yr~^ = 18.9. (Actually, only the quantity M/Mfj matters for the shape of the energy 
distribution. This, along with a*, determines the run of Tgff with radius, and gz is determined 
by T(,ff in our models. The overall luminosity can be scaled to fit the observations by scaling M. 
The same is true for the "toy" models of Section 3.) The resulting energy distribution is shown in 
Figure 2 for observer viewing angle fio = cos(0o) = 0.35. (Larger iJ,o give weaker polarization; the 
same model for /Xq = 0.1 is shown in Figure 3.) The energy distribution is in fairly good agreement 
with the observations, as found by BA and Webb et al. (1993). However, the polarization rise is 
more blueshifted and more gradual than observed. This differs from the results of BA, who did not 
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include the relativistic transfer function (but who warned that it could modify their conclusions). 

Note that for fio = 0.35, more edge on than BA's model, the polarization still reaches the observed 
~ 5 percent, but at a wavelength shorter than observed. The transfer function diminishes the peak 
polarization and adds its own blueshift to that already present in the stellar atmosphere spectra, 
causing the total blueshift to exceed that observed. 

Figure 3 shows the spectrum of the same disk, viewed at a more edge on angle (//q = 0.1) and 
scaled to the observed energy distribution of PG 1630+377. The energy distribution fits fairly well. 
However, even for this nearly edge on viewing angle, the polarization does not approach the value 
~ 20 percent observed, and the polarization rise is again too much blueshifted in the model. 

Another difficulty with the BA model may be the small ionizing luminosity from such a cool 
disk. The expected H/3 equivalent width for our model is only ~ 25 A for /Xq = 0.35 if the line 
emiting gas absorbs all the ionizing photons. Observed equivalents widths are somewhat larger, 
and the covering factor for the emitting gas may be small. Moreover, the disk's He I and He II 
ionizing luminosities are negligible. However, a "nonthermal" ionizing source may also be present. 

We conclude that the stellar atmosphere phenomenon suggested by BA to explain the observed 
polarization rises fails on account of the relativistic transfer function. The actual stellar atmosphere 
effect is nevertheless present, and it may contribute at some level to the observed spectrum of AGN 
disks. 



3. A Polarization Jump "At the Edge" 

The preceeding results show that the relativistic transfer function causes a substantial blueshift 
of any polarization rise in the disk's local emission. This raises the question, what will be observed 

if the disk gas emits a continuum with an abrupt increase in polarization right at the Lyman edge 
in the rest frame of the gas? We have calculated a simple model in which the disk surface, at 
any radius, emits a continuum with a sharp Lyman edge in absorption in the total flux together 
with a sharp jump in polarization at the Lyman edge. The emitted astrophysical flux was given by 
Fi, = 4Hi, = B^iT^ff) for v < uh and Fn = Bi,{0.84lTf,f f) for u > uh, simulating an atmosphere 
whose brightness temperature drops to the boundary temperature in the Lyman continuum. The 
polarization was assumed to rise from effectively zero at ly < i'h to an arbitrary multiple of PesiP') 
at u > uh, where pes refers to the polarization of electron-scattering atmosphere (Chandrasekhar 
1960). To this end, we define ap by p{^) = OpPesifJ-)- 

Figure 4 shows the resulting fit to PG 1630+377 for several values of the black hole angular 
momentum and /i„ = 0.1. For each value of a*, the values of M9 and Mq were adjusted to fit the 
energy distribution, and Op was adjusted to fit the maximum polarization observed. Figure 4 shows 
that the observed polarization rise remains quite sharp and is displaced to the blue by an amount 
that depends on a*. For larger a*, the disk's inner boundary at the marginally stable radius shifts 
inward, from r* = 6 for oq = to r* = 1.22 for a* = 0.9982 (Page and Thorne 1974). For larger 
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a*, the Lyman continuum comes from smaller radii where the relativistic effects are stronger. The 

wavelength and shape of the polarization rise agree with the observations for a* = 0.5. Evidently, 
this simple model gives a good fit to the observed energy distribution and to the shape of the rise in 
polarization. 

Why docs the polarization remain low for a considerable wavelength interval to the blue of 
the Lyman edge, and then rise abruptly? The relativistic transfer function causes the approaching, 
blueshifted part of the disk to appear brighter than the receding portion. In addition, light from 
the sideways moving portion of the disk behind the black hole is bent as it passses near the hole on 
its way to the observer. Even for a fairly edge-on viewing angle, this light leaves the atmosphere 
fairly close to the normal, so that its intensity is enhanced by limb darkening, the projected solid 
angle subtended by unit surface area is increased, and its polarization leaving the disk surface is low 
(cf. Agol 1997a). This light also suffers depolarization from relativistic rotation of the polarization 
plane, which varies rapidly with position on the part of the disk behind the black hole (Connors, 
Piran, and Stark 1980; Agol 1997a). The light from the part of the disk in front of the hole is more 
strongly polarized but less intense. Thus, the observed light at small Doppler shifts is dominated 
by the back part of the disk and is weakly polarized. At larger blueshifts, the light is dominated by 
the approaching part of the disk. This light is strongly polarized because it suffers less gravitational 
deflection and therefore emerges from the atmosphere at a large angle to the normal. These points 
are illustrated in Figure 5, which shows the received flux and polarization for an annulus that emits 
at a single radius, r* = 10.52, and wavelength, A = 912 A. The received emission from the near 
and far sides of the ring arc shown separately. The flux is much stronger at the blue end than at 
the red end of the observed range of wavelength shifts because of Doppler boosting. Intense but 
weakly polarized flux is received from the back side of the disk with relatively little shift from the 
emitted wavelength. Thus, the Lyman edge in the total flux will ramp up or down progessively 
from the red to the blue side of the nominal wavelength, whereas the polarization will rise only at 
shorter wavelengths. 

The blueshift of the polarization rise in the toy model depends on fig as well as a*. Figure 
6 shows results for a maximally rotating hole with a* = 0.9982. Again, M and Mh are adjusted 
in each case to flt the energy distribution. The blueshift of the polarization rise increases with 
decreasing fio , and by eye we estimate a best flt /io ~ 0.33. However, in this case, large values 
of the intrinsic polarization, Cp ~ 19, are required to match the observed degree of polarization 
of PG 1630+377. This corresponds to p(//) = 76 percent at n = 0.33 and p{n) > 100 percent at 
/J, < 0.2. The unknown mechanism for producing the postulated Lyman continuum polarization 
could have a limb darkening law quite different from pes(A*)) but polarization of roughly 76 percent 
at /X « 0.33 will likely be required in any case and may be difficult to achieve. If one rejects such a 
large intrinsic polarization, then a* must be less than 0.9982 in the context of this toy model. 

The cases /Xo = 0.25 and 0.50 are shown in Figures 7 and 8. For Ho = 0.25, the blueshift is 

best fit for a* ~ 0.7. For /io = 0.5, the blueshift depends little on a* (even for a* = 0.9982; see 
Figure 6), and it is too small to flt PG 1630+377. Figures 4 and 6-8 show that the wavelength 
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of the polarization rise provides a constraint in the (a*, /x) plane. Interpolating, we find = 
(0.10, 0.25, 0.30) for a* = (0.50, 0.75, 0.90). From the fit to the energy distribution, we find Mq 
= (27, 25, 20), Mg = (5.0, 9.0, 15), and L/Le = (0.22, 0.14, 0.099), respectively. Note that M 
decreases with increasing a* because of increasing efficiency of energy production. Also, M/M^j 
decreases with increasing a* because of stronger Doppler boosting (requiring lower Tg// ) and a 
smaller radiating area in units of {GM/c^)^, which gives higher T^ff for a given M /Mfj. All of the 
values are comfortably below the Eddington limit, consistent with a thin disk. 

This toy model was also used to fit PG 1222+228. The blueshift of the observed polarization 
rise is somewhat less than for PG 1630+377 (Impey et al. 1995; Koratkar et al. 1995), and the 
rise possibly is more abrupt, but the uncertainties are large. Figure 9 shows results for /Xo = 0.25. 
The fit for a* = 0.5 is fairly good, requiring ap = 1.7, Mq = 43, Mg = 6.2, and L/Le = 0.26. 
Reference to Figure 7 suggests that /Xq = 0.5 would fit PG 1222+228 fairly well for appropriate a^, 
independent of a*, although the blueshift may be slightly less than observed. 

Figures 6, 7, and 8 indicate that the wavelength at which the polarization rises to one-half of 
its full value, X1/2, will be between about 640 and 830 A, at least for inclinations jio < 0.5, which 
are most likely to give strong polarization. This is consisstent with the two objects discussed above 
and with PG 1338+416 (Koratkar et al. 1995). 

The observed degree of polarization for the toy model scales with ap, and results can be derived 
from the figures. For example, if we fix a* = 0.5 and ap = 4.4 (above discussion of PG 1630+377), 

we may ask what will be the plateau value of the polarization, Pmax, for this disk viewed at various 
angles. From Figures 4, 7, and 8, we find Pmax = (20, 13, 5.0) percent and A1/2 = (750, 770, 810) A 
for iJ,o = (0.1,0.25,0.5). 

4. Discussion 

We have shown that the relativistic transfer function precludes an explanation of the observed 
polarization rises in PG 1222+228 and PG 1630+377 in terms of the stellar atmosphere effect pro- 
posed by Blaes and Agol (1996). Relativistic effects add an additional blueshift to the polarization 
rise that makes the predicted rise occur at wavelengths shorter than observed. The model also 
cannot explain the strong degree of polarization observed in PG 1630+377. 

As a step toward a more successful model, we have shown that the observed polarization rises 

are consistent with a simple model in which the polarization rises abruptly at A912 in the rest 
frame of the orbiting gas. The observed polarization rise is blucshiftcd by an amount that depends 
on 0* and /v^ . For PG 1630+377, the observations are consistent with a* of about 0.5 or greater 
but probably less than 0.9982. 

What could be the origin of the polarized Lyman emission? In PG 1630+377, the polarized 
flux, F^p = PiiF^ rises steeply with decreasing wavelength (Koratkar et al. 1995). An increase 
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in polarized flux abruptly at the Lyman edge in the rest frame of the gas suggests free-bound 
emission from ionized hydrogen that undergoes electron scattering, either in the emitting gas or a 
separate scattering location (cf. Impey et al. 1995; Koratkar et al. 1995). In view of the evidence 
of irradiation of AGN disks by a hard incident continuum, at least in lower luminosity objects (e.g., 
Clavel et al. 1992; Mushotzky, Done, and Pounds 1993; and references therein), one naturally 
thinks of a photoionized layer on the disk surface (e.g., Sincell and Krolik 1997). A simple picture 
would be a uniform slab of ionized gas in which recombination of H I produces a fairly uniform 
emissivity in Lyman continuum photons that scatter off electrons in the slab on their way out of 
the slab. Such a surface layer, overlying a photosphere with a Lyman edge in absorption in its total 
flux, might resemble our toy model. Results by Phillips and Meszaros (1986) show polarizations 
of 15 percent or more for /x 0.1 in a scattering slab with uniformly distributed emitters, and the 
possibility of even larger polarizations for emitters concentrated near the surface. The strongest 
polarization occurs for Tg^ ~ 0.2 to 0.3. However, it is not clear from their results that strong 
enough polarization for PG 1630+377 will be possible, and the limb darkening law is very different 
from PesilJ')- 

The temperature of such a slab can be constrained by the fact that the polarized flux in PG 
1630+377 jumps by an order of magnitude across the polarization rise. If the polarized flux has the 
same energy distribution as the radiation feeding into the scattering process, then the emissivity 
of the gas must have a similar contrast. The emissivity ratio ji/(A912— )/j,y(A912+) decreases with 
increasing temperature because of free-free and Balmer continuum emission. The observed polarized 
flux requires T < 80, 000 K for PG 1630+377, based on expressions by Brown and Mathews (1970). 

A layer that is mechanically heated rather than photionized may be an alternative source of 
the polarized Lyman continuum. Yet another possible source involves the fact that for luminous 
QSOs, a standard (Shakura and Sunyaev 1973) disk may become optically thin in the inner radial 
zone. 

Further work is needed to establish whether the continuum and line emission in such a model 
is consistent with all known observations for QSOs with Lyman polarization rises. Koratkar et al. 
find a polarization rise at the expected position of Lyman a, weaker than the Lyman continuum 
polarization and not blueshifted. This may come from larger disk radii with lower orbital velocities. 
Simple estimates, based on the black body limit for the La surface brightness, suggest that the 
disk radii considered here will not be a significant source of La emission. 

The spin rates of black holes in AGN have received increasing attention in recent years. An 
interesting "spin paradigm", in which radio loud objects have rapidly rotating holes (Blandford 
1990) has been developed by Wilson and Colbert (1995). Values a* ~ 0.5 could arise in several 
ways. (1) If two black holes coallesce, the resulting hole has a* ~ 0.5 for a mass ratio > 0.2 if 
roughly half the angular momentum of the binary at the last stable circular orbit is radiated as 
gravitational waves during the coallescence (Wilson and Colbert 1995). (2) Successive accretion 
events at random inclinations could cause the hole's angular momentum to wander around a value 
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a* 0.5 if the mass accreted in one event is not too small as a fraction of the black hole mass 
(Moderski, Sikora, and Lasota 1997). (3) Angular momentum extraction by the Blandford-Znajek 
(1977) mechanism may balance angular momentum supplied by accretion to give an equilibrium 
value a^^eq- Results by Moderski and Sikora (1996) indicate a*^eg ~ 0.5 for am « 0.01, where a 
is the viscosity parameter (Shakura and Sunyaev 1973) and rh = c?M /LEdd- Our model for PG 
1630+377 has m 1, and values a 10~^ are not precluded in AGN. However, a value of a* as 
large as 0.5 might correspond to an object with substantial radio luminosity in the "spin paradigm", 
whereas PG 1630+377 and PG 1222+228 are radio quiet objects. 

The toy model discussed here essentially addresses the question, what must the polarization be 
in the rest frame of the gas, in order for the observed polarization rise to match the observations? 
The remarkably simple answer is that it should be an abrupt rise in polarization and polarized flux 
at the Lyman edge. This may provide a significant clue to the source of the polarized radiation. 
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Captions for Figures 

Fig. 1. — Emergent intensity, polarized intensity (dashed), and polarization at angle // = 0.1 from 
an LTE atmosphere with Tg// = 25,500 K and gz = l-2gEdd = 10^-^^ cm s~^. This T^ff corresponds 
to a radius r* = 9.5, where Tgff is highest, in a disk with mg = 5.3 and Mq = 18.9, used to model 
PG 1222+228. 

Fig. 2. — Accretion disk model compared with observations of PG 1222+228. Observations are 
from Impey et al. 1995 (filled circles), Bechtold et al. 1984 (open circles), Wampler & Ponz 1985 
(open square), Webb et al. 1993 (filled squares), and Stockman et al. 1984 (open triangles). Model 
disk with Mg = 5.3 and Mq = 18.9 is viewed at angle fio = 0.35. Position angle of the disk axis 
is unknown, so there is an arbitrary offset between model and observed position angles. Abscissa 
is rest wavelength. The position of the Lyman edge redshift is indicated in the top panel. Fx has 
been scaled to fit the observations, corresponding to modest adjustments in Mh and M at fixed 
M/Mfj (see text). 

Fig. 3. — Accretion disk model compared with observations of PG 1630+377. Observations are 
from Koratkar et al.(1995). Model has been scaled at constant M/M^ from that for PG 1222+377 
so as to fit the observed maximum flux, leaving Tf,ff{r^) unchanged. Viewing angle is jio = 0.1. 
Position angle of the disk axis is unknown, so there is an arbitrary offset between model and 
observed position angles. 

Fig. 4. — Toy accretion disk model compared with observations of PG 1630+377 (see text). Fluxes 

have been scaled to fit maximum observed Fx. Viewing angle is /i^ = 0.1. For each of the three 
values of black hole angular momentum, a*, the value of M/M^ has been adjusted to fit the observed 
energy distribution. The wavelength of the polarization rise is increasingly blueshifted for larger 
a*. The models have p/pes = (4.4, 4.4, 5.7) for a* = (0.0, 0.5, 0.9). 

Fig. 5. — Predicted flux, polarization, and positional angle for a narrow ring orbiting a black hole 
with a* = 0.5 at a radius = 10.52. The observer viewing angle is = 0.1. The emitted 
spectrum is a delta-function at A912, with the intensity and polarization in the rest frame of the 
gas depending on angle in the manner of a pure electron scattering atmosphere (Chandrasekhar 
1960.) The light received from the far side and the near side of the disk are shown separately as 
the solid and dashed curves, respectively. 

Fig. 6. — Toy model results for a* = 0.9982 for several viewing angles, compared with observations 
of PG 1630+377. The models have p/pes = (9.0, 14.9, 30.8) for /x^ = (0.1, 0.25, 0.5). 

Fig. 7. — Toy model results for jio = 0.25 and several values of a*, compared with observations of 
PG 1630+377. The models have p/pes = (6.4, 7.1, 9.6) for a* = (0.0, 0.5, 0.9). 

Fig. 8. — Toy model results for Ho = 0.5 and several values of a*, compared with observations of 
PG 1630+377. The models have p/pes = (14.0, 17.5, 24.0) for a* = (0.0, 0.5, 0.9). 
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Fig. 9. — Toy model results for fig = 0.25 and several values of a*, compared with observations of 
PG 1222+228. Values of Qp are one fourth as large as for the corresponding models applied to PG 
1630+377 (Figure 7). 
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